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Abstract 

We study the impact of modifying the vector sector of gravity on the CMB polar- 
ization. We employ the Einstein-aether theory as a concrete example. The Einstein- 
aether theory admits dynamical vector perturbations generated during inflation, leav- 
ing imprints on the CMB polarization. We derive the perturbation equations of the 
aether vector field in covariant formalism and compute the CMB B-mode polarization 
using the modified CAMB code. It is found that the amplitude of the B-mode signal 
from the aether field can surpass the one from the inflationary gravitational waves. 
The shape of the spectrum is clearly understood in an analytic way using the tight 
coupling approximation. 

Motivated mainly by the mystery of the dark components of the Universe, modification of gravity law 
from standard general relativity (GR) has been explored much. It is commonly made by adding an extra 
scalar degree of freedom as in f(R) theories. It is also possible to modify the spin-2 sector as in massive 
gravity and bi-gravity theories. In this paper, we are going to consider a hypothetical vector degree of 
freedom of gravity. Specifically, we shall focus on the Einstein-aether (EA) theory proposed by Jacobson 
and Mattingly [2] , in which a fixed norm vector field with a Lorentz- violating vacuum expectation value 
takes part in the gravitational interaction. 

The purpose of the present paper is to clarify the impact of the aether vector field on the CMB 
polarization. The CMB polarization arises from all the three types of cosmological perturbations, i.e., 
scalar, vector, and tensor perturbations. Among them, the vector perturbations most effectively generate 
the B-mode polarization [3] , though the effect has been less investigated because the vector mode decays 
unless sourced, e.g., by topological defects or the neutrino anisotropic stress [4]. Modifying the vector 
sector of gravity would add yet another possibility of producing vector perturbations, leaving a unique 
signature in the CMB polarization due to nontrivial dynamics of the aether field. 

The action of the EA theory is given by 



l - J d 4 x^[TZ - Cl \7 a A b \7 a A b - c 2 (\7 b A b ) 2 - c 3 \7 a A b \7 b A a 

- Ci A a A b V a A c V b A c ] + \{A b A b - 1)] + S m , (1) 



where k = 8ttG, 1Z is the Ricci scalar, A a is the aether field, and S m is the action of ordinary matter. 
Variation with respect to A a yields the equation of motion for the aether and variation with respect to 
the Lagrange multiplier A gives the fixed norm constraint, A a A a = 1. In the rest of the paper we use the 
following abbreviations: C13 = c\ + C3, C14 = c\ + C4, a = c\ + 3c2 + C3. 

To describe background cosmology and the evolution of vector perturbations, we employ the covariant 
equations obtained by the method of 3 + 1 decomposition. We begin with splitting physical quantities 
with respect to observer's 4- velocity u a . Following the usual procedure, the projection tensor is defined 
as h ab '■= g a b — u a u b . We define time derivative as T^...' := u c V c T b .. '. and covariant spatial derivative as 
D a T b ;;; := hfh b ■ ■ ■ h k c ■ ■ ■ V l X^ The energy-momentum tensor for each matter component and V ' a u b 
are expressed respectively as 

T ab = P {l)u aU b - p {l) h ab + 2q^u b ) + 7T^, V a Ufc = ^0h ab + (T ab + L0 ab - u a ii b . (2) 
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The energy-momentum tensor for the aether, T^ b , is also written in the same form as above. The 
expansion 8 may be written as 9 = 3S/S, where S is the averaged scale factor. 

At zeroth order, A a = u a , so that the energy density and the pressure of the aether are given 
respectively by p^ = a0 2 /6 and = — a(20 + 9 2 )/6. The Friedman equation is thus given by 
3H 2 = nS 2 p/(l - a/2) where p is the total energy density of ordinary matter and we have introduced 
the conformal Hubble parameter, H := S9/3. The background effect of the aether is just to rescale the 
gravitational constant k. 

Let us move on to the dynamics of vector perturbations. We choose u a to be hypersurfacc orthogonal, 
so that iib = at linear order. At this order, the aether field can be written as Ab — u& + DbV^ + Vb, 
where corresponds to a scalar perturbation which we do not consider in this paper, while Vb a vector 
perturbation that satisfies DbV h = 0. Each perturbation variable can be expanded using the transverse 
eigenfunctions Q± .: V a = T, V Q^ 1a ] = Eg«Q± o ab = £(fc/S)crQ± , and 7$ = E^Q^ where 
k is the eigenvalue. It is convenient to write the relevant equations in terms of the coefficient functions 
V, gW, .... 

From the equation of motion for the aether, we obtain the evolution equation for V: 

ci4 [V" + 2UV' + (H 2 + W) V]+a (H 2 - W) V + Cl k 2 V = -~^-k 2 a. (3) 

This shows that the fluctuation of the aether obeys the wave equation which is similar to the evolution 
equation for cosmological tensor perturbations. The crucial difference is the effective mass term which is 
dependent on the expansion rate % and the model parameters. The fluctuation of the aether leads to the 
effective heat-flux vector qi A ^ = — ci^D h [a a b + DuVm], which gives rise to the additional contribution to 
the momentum constraint equation. The momentum constraint under the influence of the aether is thus 
given by 

k 2 a = ( 2kS 2 Y g« - c 13 k 2 V ) . (4) 

l + c 13 \ ^ J 

Here, is determined through the individual matter equations. In calculating the CMB power spectrum 
we simultaneously solve the equations of motion for other ordinary matter and the multipole moment 
equations for photons and neutrinos, as well as the equations derived above. 

We now fix the initial condition for each variable at the early radiation-dominant epoch. This is done 
by a series expansion in terms of the conformal time 77, following [4], but now taking into account the 
presence of the aether. Neglecting the 0(k 2 ) terms, the perturbed equation of motion for the aether can 
be solved to give 



V = A^[l + Om, .:= - 1 + ^ 2 - 8Q/C14 , (5) 

implying that V can grow on superhorizon scales. Here, we have dropped the decaying mode solution 
which is not regular at r\ — > 0. The coefficient Ak is determined from the primordial spectrum of the 
aether fluctuation [5]. Requiring that scalar isocurvature modes do not grow, we consider the range 
< v < 1 [5]. As for the other variables, the appropriate early time solutions are found to be 

-^—AkV", fl «= g <*W> =0j ~- : ■W*", (6) 



1/* + 4R* 1 + ci 3 ' ' 15(1 + v) v* + 4R* I + C13' 

where we defined R* := [(1 - a/2)/(l + c 13 )]p^ / (pM + p (»)) an d v* := (5/2)(l + v)(2 + u), with the 
superscripts 7, b, and v denoting photons, baryons, and neutrions, respectively. 

In the AE theory, primordial vector perturbations are generated quantum mechanically during infla- 
tion. Once the inflation model and the reheating history are specified, one can determine the primordial 
spectrum of the vector perturbation and hence Ak, as discussed in [5]. We separate the issue of the 
perturbation evolution during inflation from the subsequent evolution, and set simply Ak = Aok^" v ~ 3 ^ 2 
(i.e., the primordial spectrum Vy cx k 71v ), where Ao is a constant. 

We have completed the numerical calculation using all the ingredients derived above and the CAMB 
code [6] modified so as to incorporate the presence of the aether. An example of our numerical results 
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Figure 1: CMB B-mode polarization and temperature anisotropy power spectra in the EA theory. For 
comparison, those from the tensor perturbation in standard general relativity are also plotted in the 
case of the tensor-to-scalar ratio r = 0.1. In this figure, c\ = — 0.2,ci 3 = — 0.3,ci4 = —a = —0.2, and 
dimensionless primordial power spectra are Vy oc fc"" and Vt oc k°. 



is presented in the Fig. 1. For comparison, we show contributions from inflationary gravitational waves 
in GR, assuming that the tensor-to-scalar ratio is given by r = 0.1. The amplitude Aq is adjusted so 
that the low-^ TT spectrum from the vector perturbation has the same magnitude as this primordial 
tensor contribution. We see in this case that the BB spectrum in the EA theory is larger than that 
from primordial tensor modes at £ > 100, and hence the B-mode is potentially detectable in future CMB 
observations aiming to detect r = O(0. 1) — O(0. 01). 

Let us try to understand the shape of the B-mode angular power spectrum Cf B in the EA theory in 
an analytic way. We start with the integral solution for the moment Bj of the B-mode polarization in 
the covariant formalism: 



e-i 
'I+i 



drjfe T ^/ e (x)C, x = k(rj - rj), 



(7) 



with ^i(x) = £je(x)/x and C = (3/4)72 — (9/2)^2- Here, ji{x) is a spherical Besscl function, r is the 
optical depth, Ig is the angular moment of the fractional photon density distribution, and Ef is the 
moment of the E-modc polarization. Using the tight coupling approximation, we obtain 

Ak Ci 3 



r 4fc 

C — a 

15f 



-V 



15f 1 + Ci 3 

It turns out that ignoring in Eq. (4) is a good approximation. We thus arrive at [5] 
V" + 2UV + c*k 2 V -] 



1 



— \H 2 

C\A, 



1 



C14 



y~o, ci = — 



Cl4 



1 



1.3 



2ci(l + ci 3 ) 



(8) 



(9) 



On supcrhorizon scales we find V oc S v in the radiation-dominant stage, as already derived, and V oc 
S v m/2 w ith v m = (—3 -f- yj\ — 24a/ci4)/2 in the matter-dominant stage. On subhorizon scale, V simply 
decays similarly to tensor perturbations, V oc S^ 1 . These relations can be mapped into the wavenumber 
dependence of V at recombination (rj = r) Tec ) as V oc Ak (k < l/c v r] Iec ), V oc k~ 2 ~ Um Ak {^/c v ri rcc < k < 
1/cu^cq), and V oc k~ x ~ h 'Ak (l/c t ,^ C q < k), where ?7eq refers to the radiation-matter equality time. 

The CMB B-modc power spectrum is roughly expressed as Cf B ~ J Vy(k)Bj Bj dhik. Using the 
approximation fe~ T ~ 5(i]-r] TCC ), Bj (r) ) can be written as Bj (m) — {k/T)[ciz/ (l+ci^V^cc)^ i [Hvo — 
Vrcc)}- (vo is the present time.) Since the projection factor ^e(x) has a peak at I ~ x, the angular power 
spectrum reduces approximately to 



C? B ~ [ ^ 



V 
Ak 



k=e/(m>-v r ,c) 



k 2+n ^^ e (k(rjo-r, Iec ))] 2 dlnk. 



(10) 
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Figure 2: (a) Scaling for the illustrative case with n v = 1 and a = — C14 = 0.2. The other parameters are 
given by C13 = —0.3, and C\ = —0.1; (b) Parameter dependence of the spectrum. In the two examples C14 
is different while the other parameters are fixed as n v = 1, C13 = —0.3, and c± — —0.1. The primordial 
amplitudes are arbitrary. 

For example, for n v = 1 and C14 = —a, the above integral can be evaluated to give the scaling £(£ + 
\)Cf B oc £ 3 (£ < 4ec := Oft - Vrcc)/c v r] ICC ), £~ 3 (£ Tec < t < 4 q := (?7o - ?7oq)/c„77 cq ), and i' 1 (£ eq < I), 
showing a peak at £ pea k ~ ^roc- This behavior can indeed be seen in Fig. 2(a), though the scaling at 
I > £ cq is hidden by the other effect and hence is not obvious. (Here, for comparison, the B-mode 
spectrum from tensor perturbations in standard GR are also plotted.) 

We can also gain an understanding of how the shape of the angular power spectrum depends on the 
model parameters. From Fig. 2(b) one can confirm the following three things: (i) since the angular power 
spectrum on the largest scales £ < £ rcc depends only on the primordial spectrum, the plotted examples 
show the same scaling; (ii) the peak position is inversely proportional to the sound velocity of the aether 
vector perturbation c v ; (iii) the difference of the small scale scaling arises due to the difference of the 
growth rate of V on superhorizon scales. The detailed discussion on the analytic estimate will be provided 
in [1]. 
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